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Abstract 

The kinetics of a 2A + B2 = 2AB reaction on supported metal 
catalyst with spillover effects is investigated using Dynamic Monte 
Carlo simulations. In the presented model A particles can adsorb 
reversibly on both metal clusters and the support whereas B2 par- 
ticles can adsorb dissociatively only on metal clusters. Particles A 
can diffuse from the support to metal clusters and vice versa. The 
model describes CO oxidation on Pd/Ai203 catalyst. Steady-state 
reaction rates and surface coverages are investigated as functions of 
partial pressures of reactants in the gas phase. Simulations show that 
spillover effects change the reaction rate dramatically by increasing it 
and widening the pressure window of the reaction. In systems with 
spillover the reaction driven by an adsorption channel, in addition to 
a spillover channel, is observed. The character of the dependence of 
the reaction rate on reactants' partial pressures can be controlled by 
changes of spillover and adsorption channel contributions. 
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1 Introduction 



Spillover is an effect accompanying heterogeneous surface reactions when 
inactive surface regions significantly influence the kinetics of the overall cat- 
alytic process. It is caused by the fact that those parts of the surface which 
are inactive in the surface reaction can be active for other processes occur- 
ring during the catalytic process, i.e. adsorption-desorption processes, and 
increase or decrease concentrations of either substrate or product particles on 
active parts of the surface through the surface diffusion. The spillover effect 
is observed and plays a significant role in catalytic reactions on supported 
metal catalysts 

The description of process dynamics in systems with spillover is compli- 
cated for several reasons. Firstly - surface heterogeneity. In most cases the 
heterogeneity, because of the occurrence of local effects, makes it impossible 
to use methods based on global chemical kinetics. Secondly - surface hetero- 
geneity means not only differences in energetic properties of different surface 
regions but also involves geometrical heterogeneities. There are active centers 
with different energetic properties on the surface but additionally those cen- 
ters are often arranged in different ways (for example active centers can form 
various clusters placed on the surface in various ways or active centers can be 
spread randomly all over the surface). The chemical and physical properties 
of the surface depend also on the arrangement of active centers. Thirdly - 
sizes of surface structures and length scales characteristic for spillover pro- 
cesses are tenths or hundredths of nanometers. Thus it is very complicated 
to use quantum chemical or molecular dynamics calculations to investigate 
such large systems and processes occurring in such long time scales. To 
investigate surface processes particularly in case of inhomogeneous surfaces 
and non-equilibrium conditions one can use Monte Carlo methods 0; [| (Ki- 
netic Monte Carlo or Dynamic Monte Carlo). Such simulation methods make 
it possible to investigate systems much larger then in the case of quantum 
chemical calculations and analyze processes with characteristic times much 
longer than in molecular dynamics simulations. These methods also let to 
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include the influence of local effects occuring on the surface. 

The current state of theoretical research of reactions with spillover effects 
includes both computer simulations and attempts at the analytical descrip- 
tion of such processes' kinetics. In one of the most recent theoretical publi- 
cations in this subject Zhdanov and Kasemo [4] derived analytical formulas 
to describe dynamics of a model reaction with spillover effect - oxidation 
of CO on a supported metal catalyst (for example Pd/A^Os). Oxidation 
reactions of this type are good examples of processes with spillover from a 
practical point of view - it is relatively simple to investigate them experi- 
mentally and they are important as actual catalytic processes. However, an 
analytical description of considered process requires many assumptions. 

Libuda, Hoffmann et al. 0; |f| presented experimental results and com- 
puter simulations of kinetics of CO oxidation on Pd/Al 2 03 but in those 
papers only mean-field simulation methods have been used and because of 
that methodology spatial heterogeneities have been neglected. 

Recently, Kovaliov [7] et al. investigated 2A + E>2 — > 2AB reaction with 
spillover using Monte Carlo simulations. In their work A adsorption was 
assumed as irreversible and they studied the influence of metal particle's 
morphology on the dynamics of the process. Our results are comparable 
with the ones obtained in their paper but in our work the metal particle's 
morphology is neglected because we wanted to focus on the spillover and to 
investigate its influence on the catalytic process. 

In this paper non-equilibrium Dynamic Monte Carlo simulations were 
used to investigate the CO oxidation reaction with the spillover effect on 
a supported metal catalyst. This allowed us to analyze the heterogeneous 
process kinetics for a wide range of its parameters and to include spatial 
heterogeneities. Using our simulation techniques new types of behavior oc- 
curring in this system have been found. 
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2 System description and methodology 



In this paper we considered the dynamics of the catalytic process that can 
be presented by the following set of chemical equations: 



A (g) + o <=> A[ ads) (1) 

A (g) + * <=> A (ads) (2) 

A\ads) + * ^ \ads) + O (3) 

B 2 (g) + ** -> 2£ (a(fa) (4) 

A( 0( fe) + £(a<fe) ^5(g) + * * (5) 



where: ^4( s ), ^4(ads), ^.( a d s ) ~~ substrate A particles in the gas phase, adsorbed 
on the active metal particle and adsorbed on the inactive support respec- 
tively; *, o - adsorption sites respectively on the active metal particle and 
on the support; ** - two neighboring adsorption sites on the metal particle; 
.Efyg) ~~ substrate B particle in the gas phase; -B( a d s ) - B particle adsorbed on 
the metal particle after the dissociative adsorption of B 2 particle; ABt g -\ - 
surface reaction product particle desorbed to the gas phase. 

The catalytic process takes place on the surface of the supported metal 
catalyst. This surface consists of two different regions: the support inactive 
in the surface reaction (for example, corresponding to AI2O3) and the active 
particles (for example, corresponding to nanometer Pd crystallites) - the only 
parts of the surface where the surface reaction given by Eq. El takes place. 
It was assumed that an A particle can adsorb, desorb and diffuse all over 
the surface. B 2 particles can adsorb exclusively on the active region where 
they dissociate and the product of the dissociation (B particles) stays on 
the active region. B particles can neither desorb nor diffuse to the inactive 
support. As mentioned above, this model can describe, for example, a CO 
oxidation reaction on a supported metal catalyst Pd/Al 2 3 j^J. 

If we neglect adsorption and diffusion (Eq. [TJ and El) we would get a well 
known ZGB model of CO oxidation reaction on a single facet of metal catalyst 
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Three additional differences between these two models exist. In our case 
A adsorption is reversible, diffusion of A particles on the metal particle is 
present and the surface reaction is not immediate. In the ZGB model A 
adsorption was an irreversible process, surface diffusion was forbidden and 
the reaction occurred immediate after adsorption of A particle on the site 
adjacent to an adsorbed B particle. We will compare our results obtained for 
the system without support (i.e. without spillover) with the ZGB model in 
Sec. O 

Analytical expressions describing the dynamics of the process in question 
can be derived assuming high surface reaction rates and high dominance of 
one substrate in the gas phase 4]. The latter assumption means that the 
value of the partial pressure of A substrate in the gas phase (pa) is near or 
1 (partial pressures of substrates obey Pa+Pb 2 — !)■ From analytical results 
one can conclude the process achieves a maximum rate in the medium range 
of pressures for various sets of system parameters. The maximum rate region 
is the most important from a practical point of view. However, analytical 
expressions derived for extreme values of reactants' partial pressures do not 
allow exact description of the medium pressures range. Our Dynamic Monte 
Carlo simulations allowed us to analyze the system behavior in the whole 
range of pressure values. 

In this paper the case of well-separated active particles was considered, 
i.e. diffusion spheres of active particles did not overlap. The investigation 
of process dynamics in such system can be reduced to simulations of the 
system consisting of one active particle placed on a relatively large support 
region. In our simulations the catalyst surface was represented by a two- 
dimensional square lattice with 150x150 nodes. In the middle of the lattice 
a smaller square consisting of 50x50 nodes was distinguished and this square 
corresponded to the metal particle active in the surface reaction. The rest 
of the lattice nodes were assumed to be a support inactive in the surface 
reaction. Periodic boundary conditions were used to exclude the influence of 
the edges of simulation lattice. 
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In simulations it was assumed that above the surface there is a mixture 
of gaseous substrates A and B 2 with known partial pressures of both com- 
pounds pa and ps 2 - According to the reaction mechanism we assumed that 
particles of substrate A can adsorb with given probabilities on active and 
inactive lattice sites (PA,ads,act and PA,ads,inact respectively). Particles A also 
can desorb from each site with probabilities PA,des,act or PA,des,inact because, 
according to Eq. |2] and ^ adsorption of substrate A is reversible. To mimic 
surface diffusion, adsorbed A particles could move on the lattice. This was 
done through jumps of A particle between neighboring lattice sites with given 
probabilities. The four nearest adjacent sites to the considered one were re- 
garded as neighboring (von Neumann neighborhood). The presence of two 
types of lattice sites (active - act and inactive - inact) allowed for four differ- 
ent types of jumps which could be characterized by different values of jump 
probabilities (from active to active site, active to inactive, inactive to inactive 
and inactive to active). 

It was assumed that B2 particles can adsorb only on active lattice sites 
and that adsorption is dissociative. Such an adsorption mechanism allows 
the adsorption of B 2 particles to take place only on those unoccupied active 
lattice sites with at least one unoccupied neighboring site - B 2 adsorption 
requires the presence of two unoccupied neighboring active lattice sites. Fur- 
thermore it was assumed that, due to the dissociation, B 2 adsorption is irre- 
versible because in the presented model the probability of B 2 desorption is 
smaller than the probabilities of other processes (B 2 desorption requires, in 
addition to crossing the surface bonding barrier, additionally the meeting of 
two adjacent adsorbed B particles and also the crossing of an association pro- 
cess energy barrier). Adsorbed B particles occupied only active sites because 
they could not diffuse from the active region to the support. In the presented 
model only A particles could diffuse from the support to the active region 
and vice versa - this diffusion is the spillover. In simulations, the diffusion of 
B particles adsorbed on the active region was neglected because B particles 
cannot cross the active - inactive region border and just that process (i.e. 
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spillover) governs the overall catalytic process dynamics. Diffusion processes 
within the active region do not directly influence this dynamics. 

If two adjacent active sites were occupied by A and B particles, irreversible 
surface reaction (Eq. EJ) could take place, an AB particle was formed and it 
was immediately removed from the lattice (corresponding to the immediate 
desorption of the product to the gas phase). The presence of neighboring A 
and B particles on the active region could be caused either by direct adsorp- 
tion of A and B particles on active sites or by surface diffusion of A particles 
from inactive to active sites. Therefore A particles could be transported to 
active regions by two different mechanisms: by direct adsorption (adsorption 
channel) and by surface diffusion from the inactive region (spillover channel). 

3 Simulation algorithm 

At the beginning of each simulation, values of partial pressures of A and B2 
substrates and values of processes probabilities were assumed. Each simula- 
tion consisted of repeated identical steps. Each step started with the random 
choice of one site on the simulation lattice. 

• If the chosen lattice site was unoccupied (i.e. there was no particle 
adsorbed on it) an adsorption attempt was realized; either an A or 
B2 particle was chosen with a probability equal to partial pressures of 
substrates in the gas phase. 

— If an A particle was chosen, an adsorption attempt was realized 
with the probability P A ,ads,act or P A ,adsMact (corresponding to the 
kind of the chosen lattice site). After this attempt, and indepen- 
dently of its results, the simulation step ended. 

— If a B2 particle was chosen, one additional neighboring site of 
the current one was randomly chosen because the B2 adsorption 
required two unoccupied neighboring sites. If the chosen neighbor- 
ing site was occupied, the simulation step ended. If it was empty, 
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an adsorption attempt on both the current and the neighboring 
site was realized with the probability P B ,ads,act or P B ,ads,inact^ in the 
presented model PB,ads,inact — because B2 could adsorb only on 
the active region of the surface). After the B2 adsorption attempt, 
and independently of its results, the simulation step ended. 

• If the lattice site chosen at the beginning of the simulation step was 
occupied and it contained an A particle, one of the three processes was 
chosen randomly: the diffusion, desorption or reaction. 

— In the case of diffusion, one of four neighboring sites was randomly 
chosen and if it was occupied, the simulation step ended. If the 
neighboring site was empty, an attempt of the particle jump to this 
empty site was made with the probability depending on the kind 
of the starting and the target site. After this attempt simulation 
step ended. 

— In the case of the desorption an attempt of removing the parti- 
cle out of the lattice with the probability P A ,des,act or P A ,des,inact 
(depending on the kind of the lattice site) was made. After this at- 
tempt and independently of its results, the simulation step ended. 

— In the case of reaction, the type of the current site was checked. 
If the current site was inactive one - the simulation step ended. 
If the current site was active - one of the four neighboring sites 
was randomly chosen. If the neighboring site was either empty or 
occupied by an A particle, the simulation step ended whereas if it 
was occupied by a B particle (adsorbed in an active site due to the 
fact that every B particle was adsorbed exclusively on the active 
region) an reaction attempt was realized with the probability P r . 
If the reaction attempt succeeded, both reacting particles were 
removed from the lattice. Independently of the reaction attempt 
results, the simulation step ended. 

• If the lattice site chosen at the beginning of the simulation step was 
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occupied and it contained a B particle, the simulation step ended. B 
particles, in contradistinction to particles A, could neither diffuse nor 
desorb in our model and the possibility of the surface reaction in the 
algorithm was taken into account when A particles were considered. 

In Monte Carlo simulations some details of implementation of physical 
processes in the simulation algorithm are important only if quantitative re- 
sults are required, e.g. if the rate constants are calculated based on the values 
of processes probabilities . In the presented algorithm, for example, the 

exclusion of the surface reaction when a B particle is found, causes that the 
reaction probability should be multiplied by the factor 2 if one wants to use 
it to calculate macroscopic quantities. In this paper we are interested only in 
the qualitative description of the system so absolute values of probabilities 
do not matter, only ratios of probabilities are important. 

The first stage of each simulation was an equilibration phase. Due to the 
fact that at the beginning of each simulation the number and configuration of 
particles on the lattice in general was changing, it was necessary to wait until 
the system reaches a stationary state (i.e. until, for given values of pressures 
and probabilities, numbers of adsorbed particles of each kind and the rate 
of the product creation were fluctuating around constant values). After the 
equilibration phase an averaging phase started - average values of numbers 
of adsorbed particles and the reaction rate were calculated based on a certain 
number of simulation steps in the stationary state. 

4 Simulations and results 

The presented algorithm was used to simulate the dynamics of processes in 
the model of the catalytic process with the spillover introduced in Sec. 121 
Parameters of simulations were chosen to correspond to assumptions which 
were taken during the derivation of the analytical description of the system 
by Zhdanov and Kasemo in [4]. This means mainly that the simulations were 
carried out always in the high reaction rates regime and therefore the value 
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of reaction probability was always high in comparison with probabilities of 
other processes. The simulations were carried out for certain sets of the model 
parameters to investigate the dynamics of the system for various values of 
reactants' partial pressures. 

The results presented by Zhdanov and Kasemo show that analytical 
solutions derived for extreme values of partial pressures, i.e. starting either 
from values of pa near or near 1 and used in a moderate region of pressures 
lead to two different types of behavior, especially that there exists multiple 
steady-state and hysteretic effects may occur. To verify the hypothesis about 
a hysteresis and to find multiple stable stationary states, for each set of model 
parameters, two different series of simulations were performed. The first 
series was started with an empty lattice for a small, close to zero, value of pa- 
Than, when the stable stationary-state was reached and average numbers of 
particles were constant, results were calculated, the value of pa was increased 
by a certain interval and a new simulation was started (therefore, in general, 
at the beginning of every next simulation the lattice was not empty). The 
procedure was repeated until the value of pa was near 1. The second series 
of simulations was carried on in an analogous way but the first simulation in 
these series was started for pa value near 1 and in succeeding simulations pa 
was decreasing. In the simulations the pa value was changing in 192 steps 
with a 0.05 interval, from 0.02 to 0.98 in the first type series and from 0.98 
to 0.02 in the second type. 

The simulations were carried on using both lattices consisting of active 
and inactive sites (this corresponded to a catalyst with an active metal par- 
ticle on the support) and using lattices built of active sites only (this corre- 
sponded to the system where the catalytic process occurs only on an active 
metal particle without contribution of the support). The comparison of the 
results obtained from both types of simulations allows to investigate the in- 
fluence of the spillover effects on the considered catalytic process. In each 
simulation surface coverages of an active and inactive region and the surface 
reaction rate were analyzed. The reaction rate was calculated as the number 
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of product particles created in a Monte Carlo step (MC step). The MC step 
consisted of simulation steps described in Sec. H3 and the number of those 
steps in the MC step was equal to the number of sites in the simulation 
lattice. 

The number of steps necessary for equilibration was determined in a series 
of introductory simulations and it was equal to 2 x 10 4 MC steps. However, 
for each first simulation in each series this number had to be increased to 
10 5 because these simulations were started with an empty simulation lattice 
and in this case the equilibration took more time. 1 x 10 4 MC steps after the 
equilibration were taken to average simulation results. Based on introductory 
simulations, an adsorption probability of particles on the active region was 
determined as PB,ads,act = 0.005. Probability of diffusion of adsorbed A par- 
ticles was taken as 1.0 independently on the kind of the starting and target 
lattice site (an isotropic lateral diffusion). The probability of surface reaction 
was equal to 0.99. Probabilities of desorption of A particles from both active 
and inactive sites were determined as 0.5 and 0.05 respectively. The choosing 
of values of the two remaining probabilities in the model (probabilities of A 
adsorption on active and inactive sites: PA, a ds,act and PA,ads,inact) allows to 
control the spillover process and to investigate its influence on the overall 
catalytic process. Adsorption and desorption probabilities in each simula- 
tion were fixed, i.e. they did not depended on the local configuration on the 
lattice, therefore lateral adsorbate-adsorbate interactions changing adsorp- 
tion/desorption energy barriers were neglected. As mentioned in Sec. in a 
case of a qualitative analysis of simulations results, only ratios of probabilities 
and not their absolute values influence the dynamics of processes. Therefore, 
in the presented simulations, values of those probabilities were chosen to 
qualitatively correspond to system properties, i.e. to get the system in high 
reaction rates regime with the dynamics controlled by the surface diffusion. 
Additionally, values of simulation parameters were optimized to achieve, in 
sequence, the best agreement with analytical results, the smallest calculation 
effort and the smallest statistical errors. 
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4.1 Reaction on the active metal particle 

In this section results of simulations with lattices consisting only of active 
sites are presented. The results correspond to the case of reaction occur- 
ring only on an active metal particle without influence of inactive support 
and therefore without the spillover effect. Probability of A adsorption was 
assumed as: PA,ads,act = 0.001. The values of remaining probabilities (as 
described in the previous section) were: P A ,des,act = 0.5, PB,ads,act = 0.005. 

In the series with increasing value of pa surface reaction does not occur 
because the surface is completely poisoned by adsorbed B particles, i.e. each 
simulation which starts with small value of pa reaches stationary state with 
the surface poisoned by B particles. To investigate the mechanism of the 
poisoning we performed additional simulations for pa = 0. In the absence of 
A particles (j>a = 0), B 2 adsorption is the only possible process in the system. 
This adsorption is irreversible, however, it cannot lead to 100% coverage of 
B particles because B 2 particles adsorb on the surface with randomly cho- 
sen orientations and therefore some lattice sites remain unoccupied. As the 
adsorption of B 2 cannot occur on single-unoccupied sites, B coverage is less 
than 100%. If surface diffusion of adsorbed B particles would be present, B 
coverage would be monolayer but in our model these particles were immo- 
bile. If pa is slightly higher than 0, an additional adsorption process occurs - 
particles A can adsorb on single-unoccupied lattice sites. Each A adsorption 
act enables surface reaction which in turn can produce a pair of neighboring 
unoccupied sites. Therefore, because B 2 is still dominant in the gas phase, 
these unoccupied sites can be attacked by B 2 and the coverage of B increases. 
The presence of certain amount of A particles is necessary for B particles to 
poison the surface. Successive increasing of pa in subsequent simulations of 
the series cannot decrease the B coverage. Adsorption of B is irreversible and 
the only way to decrease the number of adsorbed B particles is the surface 
reaction with A, however, these particles are not only absent on the surface 
but also cannot be adsorbed as well due to the absence of unoccupied sites. 
In result, in each simulation during the series with increasing pa the surface 
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is poisoned by B and surface reaction cannot occur. 

Fig. ^presents the surface reaction rate (measured as a number of prod- 
uct particles created in one MC step per one lattice site) vs. partial pressure 
of A substrate for the series with decreasing value of pa- In this series re- 
action occurs above pa — 0.82, its rate reaches maximum about pa = 0.93 
and instantly decreases near pa = 1- In the range of pressures below 0.82 
reaction does not occur. This behavior can be explained by the analysis of 
surface coverages. Fig. |2] presents surface coverages of A and B particles 
vs. partial pressure of substrate A. Even though near pa — 1 adsorption of 
A dominates, the surface is not poisoned by A. The reason of this behavior 
is the reversibility of adsorption. Particles A can both adsorb and desorb, 
therefore, the steady-state coverage is the effect of competition of these pro- 
cesses. For reasons explained later in this work, the choosing of adsorption 
and desorption probabilities (PA,ads,act = 0.001 and PA4es,act = 0.5) was done 
in such a way, to get a system with relatively small steady-state A coverages. 
Accordingly, the surface cannot be poisoned by A due to the adsorption re- 
versibility and the maximum A coverage is relatively low (about 0.004) - even 
at pa values near 1 - due to the ratio of adsorption/desorption probabilities. 

Adsorption of B2 particles is irreversible, therefore below pa = 0.82 the 
surface is poisoned by B. The mechanism leading to 100% coverage is the 
same as in the simulations starting from pa near 0. But here the reason of 
effectively low A coverage just before the poisoning is the consumption of 
A adsorbate by surface reaction. Inside the reaction window B coverage is 
less than monolayer due to surface reaction and domination of A in the gas 
phase. 

Simulations of reaction on the active metal particle can be compared with 
ZGB model. Both models show B poisoning region for pa values below the 
reaction window, however, in our case A poisoning above the reaction window 
does not occur due to the A adsorption reversibility. The reaction window in 
our model is shifted toward high values of pa with respect to ZGB because 
the steady-state A coverage is lower than in the ZGB model. 
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The presented results show that for the regions of pa above 0.82 at least 
two stable stationary states exist in the system: the state in which the re- 
action occurs (reached in the 'dec' series) and the state with the surface 
poisoned by B (reached in the 'inc' series). If the simulations start with an 
empty surface, the system would always reach the stationary state with lower 
value of the surface B coverage because the final stationary state depends on 
initial conditions of the system. In the case of experimental investigations of 
such systems the situation would be the same, if one started an experiment 
with the surface not covered by an adsorbate or covered by only one reactant 
(see P; 0]), one would reach finally only one stable stationary state. In our 
simulations, starting from both low and high values of pa caused that two 
stable stationary states could be reached. 

4.2 Reaction on the active metal particle with the sup- 
port and spillover effect 

In this section results of simulations with lattices consisting of both active 
and inactive sites are presented. Those results correspond to the case of 
the reaction occurring on a supported metal catalyst where a metal particle 
is placed on the inactive support and spillover effects occur. Probabilities 
of adsorption were assumed as following: P A ,ad s , ac t = 0.001, P At ads,inact = 
0.8, P A ,des,act = 0.5, P A ,des,inact = 0.05, P B ,ads,act = 0.005. Probabilities of 
processes occurring on the active region were equal to values taken in Sec. 
14.11 Fig. El presents the surface reaction rate (measured as a number of 
product particles created in one MC step per one active lattice site) vs. 
partial pressure of A substrate. Results of two series of simulations: with 
increasing (marked as 'inc') and decreasing ('dec') value of pa are shown. 
For comparison, results for the process with no spillover (taken from Fig. ^) 
are also presented. 

In distinction from the simulations with active lattice sites only, in pres- 
ence of spillover, reaction occurs in both series of simulations (with increasing 
and decreasing values of pa)- Reaction occurs in very wide range of pa val- 
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ues and the reaction rate is about 3 times higher with respect to the case of 
active sites only. Additionally, no differences between results of each series 
('inc' and 'dec') were observed. 

As in the previous case, one can find an explanation analyzing surface 
coverages. Fig. H] presents A coverage of the whole surface - both active and 
inactive region (measured as a number of adsorbed A particles per one lattice 
site) vs. partial pressure of A substrate. Fig. El presents the coverages of 
the active region by A and B particles (measured as a number of adsorbed 
particles per active lattice site) vs. partial pressure of A substrate. 

The curve on Fig. 0] is an isotherm of Langmuir type because about 89% 
of lattice sites are inactive and on the inactive region only unimolecular re- 
versible adsorption of A occurs. The maximum A covarage (for pa near 1) is 
lower than 100% due to the reversibility of A adsorption but is much higher 
than the maximum A coverage on the active region because of the higher ratio 

Of A adsorption/deSOrption probabilities (PA,ads,inact/ PA,des,inact > PA,ads,act/PA,des,act)- 

In the presence of the inactive support, independently on the pa value, 
the active region of the surface is not poisoned by B particles. The escape of 
B particles to the inactive support cannot be the reason of this effect because, 
according to assumptions of the model, B particles do not diffuse. Therefore 
the only reasonable explanation of this effect is the diffusion of A particles 
from the inactive support to the active region and the removing of B articles 
by surface reaction. This explanation is also supported by the fact, that the 
A coverage of active region in the current model is higher with respect to the 
model without spillover (compare Fig. EJandEJ. 

Even for low values of pa and high values of B coverage of the active 
region, the reaction occurs because the A substrate is provided from the 
support due to the surface diffusion and, at least at the border between 
the active and inactive region, the reaction always occurs. It leads to the 
decrease of B coverage on the active region and therefore A particles can 
adsorb directly on active sites from the gas phase. The effect of the lower than 
100% B coverage on supported metal catalysts was observed experimentally 
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jl^ : tl3^ . In previous Monte Carlo studies the effect was either neglected or 
included by artificial assumptions in the simulation algorithm 14]. In the 
case of spillover only one stable stationary states of surface coverages was 
found. 

The comparison between Fig. ^ and Fig. |3] allows to conclude that the 
presence of the inactive support and diffusion of A particles from and to the 
active region causes three changes in the system: a new high band with a 
maximum located about pa = 0.25 occurs, the local rate maximum in the 
range of high values of pa increases and this maximum is shifted towards 
lower pa values. We conclude that the reaction in the range of high values of 
Pa is governed by the same mechanism as in the system without the inactive 
support. However, for low values of pa the reaction is supported by the 
surface diffusion from the inactive support. In other words, in the range of 
high pa the reaction is governed by the adsorption channel, and for low pa - 
by the spillover channel. These conclusions differ from ones obtained from an 
analytical description by Zhdanov and Kasemo |4J because in our case in the 
presence of inactive support a band that comes from the adsorption channel 
still exists for high values of pa- The system with the inactive support and 
the spillover keeps properties of the system without support and shows those 
properties in the case of high values of the A substrate pressure. Similar 
results were obtained by Kovalyov et al. in |7]. Our model also does not 
show histeretic behavior in the region of moderate pa values. 

We conclude, that the increasing of the local maximum appearing in the 
range of high pressures of A with the respect to the model without spillover, 
i.e. the increasing of the adsorption channel efficiency, is caused by the ap- 
pearing of empty sites in the active region due to the additional removing of 
adsorbed particles by the surface reaction driven by spillover. The additional 
empty sites increase the rate of A adsorption on the active region and there- 
fore strengthen the adsorption channel. The same mechanism is responsible 
for the shifting of the maximum toward lower pa values. The strengthening 
of the adsorption channel manifests itself not only at the position of the max- 
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imum but increases the rate of A adsorption in the whole range of pressures. 
The shifting is caused by the asymmetry of the strengthening, the reaction 
rate in the region governed by the adsorption channel is increased stronger 
for lower values of pa- 

To investigate more precisely the influence of spillover channel on the be- 
havior of both bands occurring on the rate vs. pa plot, a series of simulations 
analogous to previous one but for different values of A adsorption probabil- 
ity on the inactive region were carried on. The increasing of this probability 
leads to the increasing of the steady-state coverage of inactive region by A 
particles and increases the flux of A particles from the support onto the ac- 
tive region. Fig. |H1 presents the maximum reaction rate and the pa pressure 
corresponding to the maximum obtained in series of simulations for one set 
of model parameters (i.e. the maximum rate obtained on the plots analogous 
to Fig. El) vs. probability of adsorption of A particles on the support (all 
model parameters for simulations presented on this plot except PA,ads,inact 
were fixed). The points for PA,ads,inact = are taken from simulations of the 
system without support i.e. for the system without spillover effects. 

Fig. |U1 shows that the increasing influence of spillover channel leads to 
the increasing of the maximum reaction rate whereas the pa pressure corre- 
sponding to that maximum, is decreasing. When the spillover becomes more 
important, the reaction rate increases in the range of low pa pressures. In the 
case of weak spillover, the system acts as a system consisting of active centers 
without a support. This effect is shown on Fig. which presents the reaction 
rate vs. pa for different values of probability of A particles adsorption on the 
support. 

4.3 Spillover effects with the strong adsorption chan- 
nel 

Among various models considered in the present work (with different values 
of all parameters) it is worth to present the case where an adsorption of A 
substrate on the active region is much faster than in simulations presented 
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in the previous section. Fig. |H] presents the reaction rate vs. pressure of A 
substrate for the model with the probabilities: PA,ads,act = 0.01, PA,ads,inact = 

0.8, PA,des,act = 0.5, PA,des,inact = 0.05, Ps^ads^act = 0.005. Here, PA,ads,act 

probability is 10 times larger than in the simulations presented in previous 
sections. Reaction rate for both the system consisting of active sites only 
(i.e. without spillover effects) and for system consisting of the active region 
and the inactive support (with spillover) are presented. 

In this case, when fast adsorption on inactive region is present, the re- 
action rate is high even in the system consisting only of active sites whit 
respect to simulations presented in Sec. 14.11 The range of pa values for 
which the reaction occurs is much wider than in the cases presented above 
and the rate maximum is shifted toward lower pa values (compare with Fig. 
HJ. The presence of the inactive support slightly increases the maximum re- 
action rate and shifts the maximum left. Additionally, spillover broadens the 
range of active pressures and a new band occurs below pa = 0.3. But since 
reaction on the isolated active region is very fast, the height of the band 
located above pa = 0.3 is still bigger than the height of the band coming 
from the spillover channel. This effect is different than in the cases presented 
in previous sections because now, in contradiction to Fig. even for rela- 
tively low values of the ratio PA,ads,inact/PA,ads,act (this ratio corresponds to 
the strength of spillover channel in comparison with the adsorption channel 
and on Fig. |H] is equal 80) the reaction rate in the range of low A pressures 
increases and the reaction occurs always in the whole range of pa values. In 
previous simulations the reaction occurred for the whole range of pa values 
only if P A ,ads,inact > 0.4 (then P A ,ads,inact/ ' PA,ads,act > 400), i.e. for very strong 
spillover channel. 

Results for our model with strong adsorption channel may be compared 
with results presented by Kovalyov et al. in |2j ■ In that work the presence 
of surface diffusion from and to the inactive support causes slight increasing 
of the maximum reaction rate and a new, but not well resolved, band occurs 
for low pa values. Those results are comparable with ours only qualitatively 
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due to the differences in the values of probabilities chosen in both models. 

It is noteworthy that for certain values of A adsorption probabilities on 
both the active and inactive surface region one can get the system in which 
the reaction rate will be almost independent of reactants' partial pressures. 
Such an experimental system with CO oxidation on Pt/Sn02 catalyst was 
reported, for example, by Grass and Linz but in that work the system be- 
havior was explained using a different model (a reaction mechanism including 
an oxygen migration was proposed). 

5 Conclusions 

Supported metal catalysts show both energetic and spatial heterogeneities. 
Dynamic Monte Carlo simulations can treat those inhomogeneities explicitly. 
In this article results of Dynamic Monte Carlo simulations of the model 
reaction of CO oxidation on the supported metal catalyst are presented. 
In this reaction the influence of oxide support is significant because of the 
presence of spillover effects and therefore the importance of spillover was 
investigated. We simulated catalytic systems both with and without spillover 
effects and in the former case the importance of spillover was controlled by 
adjusting the spillover contribution. 

1. In systems without spillover: 

(a) The poisoning by A particles does not occur due to the A adsorp- 
tion reversibility. 

(b) The catalyst surface cannot be poisoned by the B reactant if re- 
actant A is not present because of the two-center dissociative ad- 
sorption of B. Only if reactant A is present and its partial pressure 
is low, B particles can poison the catalyst. 

(c) Two stable stationary states of B coverage were found 

(d) The reaction rate achieves maximum for relatively high values of 

PA- 
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2. In systems with spillover: 

(a) The surface cannot be poisoned by B particles for any value of pa- 
If pa = the effect is caused by two-center adsorption whereas for 
Pa > the reason is the constant presence of A particles on the 
active region due to diffusion from the support. We have obtained 
this effect without additional assumptions in the algorithm as in 
the previous studies. 

(b) The reaction is possible for every value of pa (except and 1) 
because the surface cannot be poisoned. 

(c) Depending on the value of pa, reaction occurs in two regimes: 
driven by the spillover channel and driven by the direct adsorption 
of A on the active region of the surface. 

(d) No hysteretic behavior have been found. 

(e) Relative values of reaction rate in the case of reaction driven by 
the spillover channel and reaction driven by the direct adsorption 
of A on the active region depend on relative values of probabili- 
ties of adsorption and desorption of A on the support and on the 
active region of the catalyst surface. For certain values of those 
probabilities one can get a system in which the reaction rate will 
be almost independent of partial pressures of reactants. 

Some of our results differ qualitatively from the ones obtained by Zhdanov 
and Kasemo in jj] for an analytical model: we have not found hysteresis in the 
system with spillover and we have found two reaction regimes i.e. driven by 
spillover and driven by adsorption because we were able to use one model for 
the whole range of values of reactants' partial pressures. Analytical solutions 
require additional assumptions about partial pressures and surface geometry, 
whereas Dynamic Monte Carlo simulations can be used for wide variety of 
model parameters and can treat surface heterogeneity explicitly. 

Our results are in good qualitative agreement with the Kovalyov and co- 
workers' results Q. Results for their model with spillover in the case of flat 
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and uniform active particle are comparable with our simulations for strong 
adsorption channel. However, our data for a strong spillover channel are 
qualitatively new. 

Dynamic Monte Carlo simulations are a useful tool to study processes 
occurring on heterogeneous surfaces, but because of the lack of experimental 
data which are needed to calculate model parameters, in most cases their 
results should be treated only qualitatively. Therefore results presented in 
this work show qualitatively what kind of effects can occur in systems with 
spillover and which parameters influence them. 
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Figure 1: Surface reaction rate vs. partial pressure of A in the gas phase for 
simulation series with the decreasing ('dec') value of pa in the case of the 
active metal particle without the support. 
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Figure 3: Surface reaction rate vs. partial pressure of A in the gas phase for 
simulation series with the increasing ('inc') and decreasing ('dec') value of pa 
in the case of the active metal particle with the support and spillover effects. 
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